CHARACTERIZATION OF RADIALLY SYMMETRIC FINITE TIME 
BLOWUP IN MULTIDIMENSIONAL AGGREGATION EQUATIONS* 
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Abstract. This paper studies the transport of a mass fj, in R d , d > 2, by a flow field v = ~\7K*fi. 
We focus on kernels K = for 2 — d < a < 2 for which the smooth densities are known to 

develop singularities in finite time. For this range we prove the existence for all time of radially 
symmetric measure solutions that are monotone decreasing as a function of the radius, thus allowing 
for continuation of the solution past the blowup time. The monotone constraint on the data is 
consistent with the typical blowup profiles observed in recent numerical studies of these singularities. 
We prove monotonicity is preserved for all time, even after blowup, in contrast to the case a > 2 
where radially symmetric solutions are known to lose monotonicity. In the case of the Newtonian 
potential (a = 2 — d), under the assumption of radial symmetry the equation can be transformed into 
the inviscid Burgers equation on a half line. This enables us to prove preservation of monotonicity 
using the classical theory of conservation laws. In the case 2 — d < a < 2 and at the critical exponent 
p we exhibit initial data in L p for which the solution immediately develops a Dirac mass singularity. 
This extends recent work on the local ill-posedness of solutions at the critical exponent. 

1. Introduction. 

This manuscript considers the problem of dynamic nonlocal aggregation equations 
of the form 

^ - div(pVK * p) = (1.1) 
in R d for d > 2. This problem has been a very active area of research in the literature 

iiiiaiaiiniiiiEaiiiiniiiiir^ 

[SniinilSllSlinilSllMllZniin]- These models arise in a number of applications 
including aggregation in materials science (4TJ [42j EH [65], cooperative control [39] . 
granular flow [25] [26] [71] , biological swarming models [62] [69] [70] , evolution of vortex 
densities in superconductors [36] El [2j [35] [58] and bacterial chemotaxis [21] 06] [TB] [T7] . 
A body of recent work has focused on the problem of finite time singularities and local 
vs global well-posedness in multiple space dimension for both the inviscid case (|1.1[) 
d EH [H H QJ HD 133 EH EE] and the cases with various kinds of diffusion 
[5] [TB] [13] [50] [51] . The highly studied Keller- Segel problem typically has a Newtonian 
potential and linear diffusion. For the pure transport problem (jl.ip . of particular 
interest is the transition from smooth solutions to weak and measure solutions with 
mass concentration. This paper presents a general framework for radially symmetric 
solutions that blowup in finite time in which the initial data decreases monotonically 
from the origin. This paper differentiates itself from the previous work in that it 
considers continuation of the solution as a measure past the initial singularity for 
the range of 1 > a > 2 — d. Prior work on measure solutions considers the case 
a > 1 in general dimension [24j and the Newtonian case with a defect measure in 
two dimensions [65] . The monotone constraint is typical of the local structure of the 
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solution near the blowup time as shown in numerical simulations. This paper presents 
rigorous theory for solutions with such structure, showing that the nonlocal evolution 
preserves this global structure. Much work has been done on related problems in fluid 
dynamics to understand measure solutions of active scalar equations [30l [31] [60l [65l 
[77] . The results of this paper may provide further insight for those problems. 

In the case of power-law kernels e.g. K{x) = \x\ a /a, it is already known that the 
critical power is a = 2. For a < 2 finite time singularities always arise and for a > 2 
solutions stay smooth for all time if the initial data is smooth [9l [TUJ [24] . In the case 
of finite time blowup, i.e. a < 2, it has been observed numerically (44l 05] 03] that, 
starting with a smooth radially symmetric initial data which is monotone decreasing 
as a function of the radius, the solution evolves so that the monotonicity is preserved 
and at some finite time develops a power-law singularity at the origin. The power is 
sufficiently singular that, based on the results in this paper, the solution produces an 
instantaneous mass concentration after the initial singularity. To make this rigorous 
we must develop a general theory for such singular solutions. This paper addresses 
this particular class of measure solutions, namely those with a radially symmetric 
decreasing profile and possibly a Dirac mass at the origin. We prove that this structure 
is maintained for all time when a < 2. We conjecture that this class of radially 
symmetric decreasing solutions including a Dirac mass at the origin describes well 
the local behavior at the blow up time of general non radially symmetric solutions. 
We also note that for a > 2 it has been observed in numerical simulations 05j 03] 
that monotone decreasing structures are not preserved: indeed there is an attracting 
solution of the form of a collapsing delta-ring which causes mass to collect on the ring 
during the collapse thereby destroying any initial monotone property of the solution 
05103]. 

We note that our work fits nicely between the general measure theory in EI] 
(for a > 1) and previous works which consider more regular classes of weak solutions 
including L°° [8] [9] and LP [12]. For L°° and LP solutions we typically have only 
local well-posedness whereas the measure solutions have global well-posedness. Our 
work extends the global existence results in [24] to the case of more singular kernels 
with power a > 2 — d, for the special case of monotone decreasing radially symmetric 
measure solutions. This includes that of the Newtonian potential, which is discussed 
separately in the next paragraph. Uniqueness of solutions for 2 — rf<a<lis still 
an open problem. 

In two space dimensions, when K(x) = log(|x|) (i.e. K is the Newtonian poten- 
tial), the aggregation equation arises as a model for the evolution of vortex densities 
in superconductors [36] |66] E3 023 El El [56l [35l [58] , and also in models for adhesion 
dynamics |64 | I65j. In these models singularities are known to appear in finite time, 
and the question of interest is how to continue the solution after the initial formation 
of singularities. Since these singularities are expected to be Dirac masses one has to 
consider measure solutions. Unfortunately, due to the very singular behavior of the 
Newtonian potential at the origin, most of the results to date concern the existence 
of measure solutions which contain an error term (a defect measure) compared to the 
original equation [35l[64][2]. Also uniqueness is lacking in these works. 

In this paper we consider the Newtonian case in all dimensions and show that for 
general radially symmetric data there is no need to consider a defect measure because 
the symmetry allows the problem to be reduced to a form of the inviscid Burgers 
equation on the half line, for which many things arc known. In particular, the case of 
radially symmetric monotone decreasing densities maps to classical Lipschitz solutions 
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of the inviscid Burgers equation, without shocks, allowing us to prove such solutions 
exist and are unique for all time. For the non-monotone case, shocks can form, 
corresponding to mass concentrations along spherical shells, and their evolution is not 
immediately well-defined, due to a jump in the velocity field at the shell. However, 
one can use the classical weak solution theory for Burgers equation to define a jump 
condition through a weak form of the evolution equation, or through some other 
convention. If we use the classical Burgers shock solution then the solution is unique 
since it automatically satisfies the Lax entropy condition. The more singular case 
of signed measures can also be studied in this framework, in which case one must 
consider rarefaction solutions as well as shocks. 

Going back to the case K(x) = \x\ a /a, 2 — d < a < 2, recent computational 
results |44j show that the initial finite time blowup from radially symmetric data 
has a simple self-similar form in which the power-laws of the similarity solution have 
anomalous scaling but the shape of the similarity solution has a simple monotonically 
decreasing structure with powerlaw tail. The power in the tail determines the degree 
of singularity of the solution at the initial blowup time - we observe that at the initial 
blowup time the solution leaves L°° but remains in some LP spaces and does not 
concentrate mass. This result prompted a careful study of the well-posedness of the 
equation in LP spaces |12j . In that paper it was proved that for a given interaction 
kernel, there exists a critical LP space such that the problem is locally well-posed for 
p > p c . Moreover it was proved in |12j that the power p c is sharp for K = \x\, i.e. 
the problem is locally ill-posed for p < p c . Some of these results, in particular the 
critical LP space, have been extended to general power- law kernels in [33]. In the 
present work we examine the mechanism by which initial data in the critical space 
LP C leave instantaneously this space. Taking advantage of our existence theory for 
radially symmetric decreasing measure solutions when 2 — d < a < 2, we exhibit a 
large class of radially symmetric decreasing initial data in LP" for which a Dirac mass 
forms instantaneously in the solution. This is a natural extension of the results in 
[H and [33]. 

This paper is organized as follows: below we review the mathematical notation 
and basic functional analysis used in this paper. Scction[2]dcvclops a general existence 
theory for radially symmetric solutions with the monotonicity constraint. Section [3] 
proves instantaneous mass concentration for the critical LP spaces. Section[?]considers 
the case of the Newtonian potential, for which we can show that radial symmetry 
results in a transformation of the nonlocal problem to the inviscid Burgers equation 
on a half line. Section [5] summarizes the results and discusses some open problems. 
In the appendix we derive some background theory of ordinary differential equations 
needed for the proofs in this paper and not derived in standard references (although 
the arguments are similar to standard methods). 

1.1. Mathematical formulations and notation. 

The aggregation equation, for smooth solutions, in Eulcrian coordinates, is 



For very singular kernels, and correspondingly singular solutions - in general measure 
solutions - it makes sense to reformulate the problem in Lagrangian coordinates and 
work mainly in this framework to develop the theory. It is easy to see, in the case of 
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— + div(pu) = 
at 

v(x, t) = —VK * p. 



(1.2) 
(1.3) 



strong solutions, that the above formulation is equivalent to 



p(t) = cr'#p™ t , (1.4) 
a is the flow map associated to the field v = —(VK * p(t))(x). (1-5) 

In other words, the mass p is transported by characteristics a that satisfy the ordinary 
differential equation 



— a[x, t) = v(<t(x, t), t) , o~(x, 0) = x. 

The map cr* :M. d — >• M. d is defined by cr' (x) = a(x, t) and <t* stands for the push 

forward of the measure Pi n u by the map cr' (see below for a precise definition of the 
push forward). We work with formulation (|1.4H1.5[) to prove existence of solutions, 
rather than (|1.2H1.3[) . We refer to this as the Lagrangian formulation of the problem. 
Note, in particular, that the flux pv in (|1.2[) may be difficult to define, for the product 
of a measure p and a velocity field that blows up precisely at the point where p 
concentrates mass. This phenomena can occur, for example, in the case of power 
law potentials potentials K(x) = \x\ a /a, a < 1, for which existence theory was not 
known prior to this work. Since we are working with a purely transport problem it is 
very natural to work in a Lagrangian framework. The radial symmetry combined with 
monotonicity provides a focusing effect in which the only mass concentration occurs 
precisely at the origin, providing a natural way to keep track of mass transport in this 
problem. We now introduce some technical notation and corresponding well-known 
functional analytic results. 

• A4(M. d ) stands for the space of Borel non- negative measure on M. d which have 
finite mass. 

• A4 R (M. d ) is the set of p G Ai(M. d ) which arc radially symmetric. 

• ■Mrd^-' 1 ) is the set of p € M(M. d ) which are radially symmetric and de- 
creasing. To be more precise, p belongs to Ai R D(^- d ) if and only if it can be 
written p = mS + g, where m £ [0, +oo), S is the Dirac delta measure at the 
origin and g is an L 1 function which is nonnegative, radially symmetric and 
monotone decreasing as a function of the radius. 

• V(R d ), V R (M. d ) andr RD (M d ) are the subset of M{R d ), M R (R d ) and M RD (R d ) 
respectively which are made of measure of mass 1. 

• T , 2 0^ d ) C P(M. d ), is the subspace of probability measure of finite second 
moment, i.e. j Rd \x\ 2 dp{x) < oo. 

• We say that a sequence (p n ) C T'(M d ) converges narrowly to p E 7 , (IR' i ), 
denoted by p n p, if 

lim / f(x)dp n (x) = / f(x)dp(x) 

for every / 6 C°(K rf ), the space of continuous and bounded real function 
defined on R d . 

• C w {[0,+oo),V{R d )) is the set of functions p : [0, +oo) -> V(R d ) which are 
narrowly continuous, i.e. p(t + h) — p(t) as h — > Vt > 0. 

• For p and v in ^(K^), Wiiji^u) stands for the Wasserstein distance with 
quadratic cost between p and v (See [65] for the definition and properties of 
the Wasserstein distance Wi{p,v))- Recall that P2(K d ), endowed with the 
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metric W2 is a complete metric space. Furthermore, 

lim W%{nn, n) = => /U n — 1 fj, as n — >• 00. 

• C([0, +oo),"P2(R d )) is the set of functions from [0, +00) to 7> 2 (R d ) which are 
continuous with respect to W2. Note that 

C([0,+oo),7? 2 (R d )) C a„([0,+oo),P(R d )). 

The space C([0, +00), ^(H^)) is endowed with the distance 

W 3 (ji,v) = sup Wa(ji{t), i/(t)). 

i>0 

• If T : R d -> R d is a Borel map, and if ^ e M(R d ), we denote by Tffp the 
push forward of fj, through T, defined by T#fi(B) = /i(T- x (B)), VB G i3(R d ). 
More generally we have 

f(T(x))dp(x) = [ f(x)d(T#»)(x) 

for every bounded Borel function / : R d — > R. 

• Both B(0, R) and £?r will be used to denote the open ball of radius R, {x G 
R d : \x\ < R}. A e , e < 1, denotes the annulus {x G R d : e < |x| < 1}. 

All the probability measures in this paper are compactly supported and therefore 
belong to the space V2 (R n ) on which the Wasserstein distance with quadratic cost is 
defined. 

1.2. Lagrangian solutions versus distributional solutions. 

Let us focus on power law kernels K — \x\ a /a for simplicity. If v is bounded 
on compact sets, which is only true for a > 1, it is then standard (see [T] or [12l 
Proposition 4.8] for example) to prove that if p and a satisfy (|1.4p and (|1.5p . then p 
is a distributional solution of the aggregation equation, i.e. 



+00 



( -£(x,t) + Vf(a?,t) -«t(s) ) d/>t(x)dt = 0, (1.6) 
: \ at / 

»,(i)=»(i,i) = -(VJ(*p(i))(i) (1.7) 

for all £ G C^°(R d x (0, +00)). On the other hand, if 2 - d < a < 1, then the velocity 
field x M> vt(x) is not bounded and it is not clear how to give a sense to (|1.6j) . Hence 
we use the Lagrangian formulation of the problem throughout most of this paper. In 
the special case of the Newtonian potential, in Section 0] we transform using mass 
variables to Burgers equation for which it again makes sense to use a distributional 
form of the problem albeit in a different coordinate system. 

2. General theory of radially symmetric decreasing solutions. 

This section is devoted to the proof of the following theorem: 

Theorem 2.1. Let VK(x) = x\x\ a ~ 2 , a G (2 - d,2). Given p mtt G T RD {R d ) 

with compact support, there exists p G C([0, +00), VRD(R d )) md a continuous map 

a : [0, +00) x R d ->■ R d satisfying 



p(t) = o-tftpinu, (2.1) 

'-(VK*p(i))(x), x^O 
0, i = 0. 



, n 1 1 nil/ \ \-(VK*P{t))(x), 1^0 ,„„, 

o~ is the flow map associated to the field v(x,t) = < (2-2) 
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Remark 1. In (|2.2|) we could have let v(x,t) = —{VK * p(t))(x) for all x with 
the understanding Vif(O) = and the convolution is taken in the sense of principle 
value. 

Theorem 12.11 is interesting for two reasons: first it provides global existence of 
radially symmetric decreasing measure solutions with potential more singular than 
the one considered previously. In [24| global existence and uniqueness of measure 
solutions is proven for a > 1; here we restrict our attention to radially symmetric 
decreasing solution but we obtain global existence for 2 — d < a < 2. Secondly this 
theorem shows that radially symmetric decreasing profiles are preserved for all time 
when 2 — d < a < 2. Monotonicity is also preserved for the Newtonian case however 
in this case the problem localizes and the simpler proof is carried out in Section 0J 

2.1. Formula for the convolution in radial coordinates and properties 
of the kernel. 

In this section we recall some known results about radially symmetric solutions 
of the aggregation equation and we prove additional results needed in the following 
subsections. 

Definition 2.2. Let n e MB.(R d ). We define fi e .M([0,+co)) to be the Borel 
measure on [0, +oo) which satisfies 

fi{I) = n({x G R d : \x\ G /}) 

for all I G B([0,+oo)). 

Remark 2. It is straightforward to check that if a sequence /J, n G Vn(M. d ) con- 
verges narrowly to ji G ■Pr(R <z ) then 

hm / f(r)dfi n (r) = f f(r)dfi(r) (2.3) 

™^°°J[0, + oo) J[0, + oo) 

for every f G C°([0, +oo)), the space of continuous and bounded real function defined 
on [0, +oo). 

Definition 2.3. Let [i, v G Vn(M. d ). We say that [i is more concentrated than 
v , and we write fi >- v, if fi = Tf^v for some Borel map T : [0, +oo) M> [0, +oo) 
satisfying T(r) < r for all r G [0, +oo). 

Remark 3. Suppose that /i, v G Vn(R d ) with v = m5 + / for some m > and 
f G L 1 (M. d ). It can then be proven that 

fi>v fi (b(0, r)) > v (B(0,r)) Vr > 0. (2.4) 

This equivalence is actually true for any fj,,v G T'R(R d ) if one uses transport plans in 
Definition \2.3\ rather than only transport maps ( see J73jj for a definition of transport 
plans). The proof is a consequence of the fact that, given any two probability mea- 
sures on [0, +oo) , an optimal transport plan with respect to the quadratic cost can be 
explicitly constructed in term of the cumulative distributions of these two probability 
measures \7S\ Theorem 2.18 page 74]- It can then be checked that if the cumulative 
distribution of fi is greater than the cumulative distribution of v, this optimal transport 
plan takes elements of mass from v and move them toward the origin. It can also be 
checked that if v = mS + / for some m > and f G -^([O, +oo)), then this optimal 
transport plan is induced by a transport map. Since the equivalence (|2.4j) is not used 
in this paper, we omit the proof. We note that other authors (see for example \72\ 
Chapter 1]) use the right hand side of (|2.4[) to define the notion of concentration. 



6 



For ae(2-d,2) define the function (ft : [0, +00) -> R by 

0(0 = — / , Cl ~T a • ei Mv), (2-5) 
w<j-i 7srf-i |ei — ry\ 

where S d ~ 1 — {x £ M. d : \x\ = 1} is the unit sphere and w<j-i its surface measure. 
The following lemma was proven in |12j for the case a = 1 and in |33j for general a. 

Lemma 2.4. Let VK{x) = x\x\ a ~ 2 , a £ (2 - d, 2). Let /i £ _M fl ,( Rd )- Then for 
any x ^ 0, we have 

{VK * /i)(a?) = M"" 1 f °° (ft (-!-) dfi(r) ^r. (2.6) 



Moreover, (ft is continuous, strictly positive, non-increasing on [0,+oo), anc 

d + a-2 



2-a 



0(0) = 1, lim (ft(r)r 

r~>oo a 

Note in particular that (ft £ Cq([0, +00)) which, in view of (|2.3j) . will be convenient 
in order to pass to the limit in expressions such as ()2.6|) . The positivity, monotonicity 
and boundedness of (ft have three important consequences that can be directly read 
from (HU). 

Corollary 2.5. Let /i, v £ -Mi?(R d ). Since (ft is strictly positive, we have 

H>v => |VX * /i| > |VJf * i/| . (2.7) 
Let /j,, v £ Vn(BL d ). Since (ft is non-increasing, we have 

pi>-v => |VX * fi\ > \VK * v\ . (2.8) 
Let fi E VR(R d ). Since < (ft < 1 we have 

\WK * n\ < l^l"" 1 . (2.9) 



In the next Lemma we prove that that (ft is C 1 for a > 3 — d, quasi-Lipschitz 
continuous for a = 3 — d and Holder continuous for 2 — d < a < 3 — d. The lack of 
smoothness for a < 3 — d is due to a singularity in the derivative r = 1. Below we 
prove sharp estimates on the regularity of (/>; later we will only use the fact that eft is 
Holder continuous in the range of a considered. 

Lemma 2.6. 

(i) If a S (3 — d, 2) then (ft £ C 1 (0, +00) and (ft 1 is bounded on (0, +00). 
(ii) If a = 3 — d, then there exists a constant C\ > such that 

[(ftin) - <ft(r 2 )\ < Ci |n -r2|(l- login -r a |) ( 2 - 10 ) 

for all n,n satisfying |n — n| < 1/2. 
fm) If a £ (2 — d, 3 — d), i/ien i/iere exists a constant C2 > suc/i 4/ia£ 

|0(n) - 0(r 2 )| < C 2 |n - r 2 \ a - (2 - d) (2.11) 



/or all r\,t2 satisfying |n — n| < 1/2. 
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Proof. In [33} Lemma 4.4] it was proven that <j> is diffcrentiablc on [0, 1) U (1 + oo) 
and that, for r ^ 1: 

where A(r,8) = (1 + r 2 - 2rcos6») 1/2 (2.13) 

and 6 Q , d = — — — > 0. (2-14) 

Wd_i(0 - 1) 

Note first that for fixed r the function 9 >-»■ A(r, 0) reaches its minimum at = 0. So 
-^('"i 6) > |r — 1| an d one can easily sec from (|2.12|) that for all the a considered, <j>'(r) 
is bounded on [0, 1/2] U [3/2, +oo). It is therefore enough to prove the statements of 
the Lemma only on the interval (1/2,3/2). 

We first prove (i). Note that for fixed 9 the function r H> A(r, 9) reaches its 
minimum at r = cos9 and therefore A(r, 9) > |sin0|. As a concequence we have the 
following estimate for the integrand in (|2.12l) : 

r(S[li9)d < „JjL- d forallrefi f) and e [0, tt]. (2.15) 



A(r,9) 4 ~ a ~ (sin9) 4 ~ a - d \2' 2 

Since the right hand side of (|2.15l) is intcgrablc if a > 3 — d we obtain (i) by the 
dominated convergence theorem. 

We now turn to the proof of (ii) and (iii) . We first derive the estimates 

\<l>'(l + h)\< ni3 C a _ d ifae(2-d,3-d) (2.16) 
\h\ 

\<f>'(l + h)\ <-Clog\h\ if a = 3-d (2.17) 

for all h £ (—1/2, 1/2), h ^ 0. The constant C > depends on a but not on h. We 
will prove (g7TS| only for h e (0, 1/2). The proof for h e (-1/2,0) is precisely the 
same. Write 



r | / " k (^ + h)(shx9) d f (l + h)(shx9) d 



and let (/) be the first integral and (II) the second one. Using the fact that A(r, 9) > 
sin 9 we find that 



fTT/2 i p 

< 3 / -ts— -d9 < 



— cx—d 



h (f0) 4 - a - d " V 

where we have used the symmetry of sin 9 around 9 = tt/2 and the fact that sin 9 > 
(2/ir)9 on the interval [0, 7r/2]. To estimate (II) we use the fact that A(l + h,9) > h: 



(II) <-^f (^9) d d9 < JjL (\ s in6fd6 

n J[0,h]U[ir-h,Tr] n JO 



h .,.. C 



This concludes the proof of (|2.16l) . The proof of (|2.17|) is similar. Let 

, n L ?■> rn r a -( 2 -V ifaG(2-d,3-d) 
w (r) = <^ Q ~ (2 - d) v ; (2.18) 

[CY(l-logr) if a = 3-d 

and note that u> is the antidcrivative of the right hand side of (|2.16[) and (|2.17|) . Note 
also that w is a nonncgative, increasing, concave function on [0, 1] which is equal to 
at r = 0. To conclude the proof of (ii) and (iii) we need to show that w is the modulus 
of continuity of on the interval (1/2,3/2), that is 

mn) - Hr 2 )\ < Lo(\ n - r 2 \) foralln,r 2 e (2.19) 

Since 0' is negative, from (|2.16|) and (|2.17[) we have that < — <f>'(l + h) < uj'(h). for 
h > 0. Integrating this inequality on [hi, h 2 ] and using the fact that w(h 2 ) — <jj{hi) < 
w(/i2 — hi) — w(0) = uj{h 2 — hi) due to the concavity of uj, we obtain that 

< 0(1 + h x ) - 0(1 + h 2 ) < oj(h 2 - hi) for all < hi < h 2 < 1/2. 

This prove that (|2.19[) holds for all ri,r 2 £ [1, 3/2). A similar proof leads to the same 
result on the interval (1/2, 1]. To obtain the result on the full interval (1/2, 3/2), let 
hi,h 2 G (0, 1/2) and write 

< <t>(l - hi) - 0(1 + h 2 ) = 0(1 - hi) - 0(1) + 0(1) - 0(1 + h 2 ) 

< ui(hi) + u(h 2 ) < 2u(hi + h 2 ). 

To obtain the last inequality we have used the fact that u is increasing. □ 
2.2. Regularity of the velocity field. 

We now study the regularity of a the velocity field associated with a radially 
symmetric decreasing measure solution of the aggregation equation. Recall that 

A e := {x eR d : e < \x\ < 1}. (2.20) 

Obviously A Q = 5(0, 1)\{0}. 

Proposition 2.7. Let p e C w ([0, +oo),V RD (^ d )) and assume that supp{p(t)) C 
5(0,1) for all t > 0. Let VK(x) = x\x\ a ~ 2 , a G (2 - d,2). Then the velocity field 
v(x, t) defined by 

^t)^-^* 9 ^ (2-21) 
v ' [0, x = v J 

satisfies: 

(PO) v(x,t) is continuous on Aq x [0, +oo) for all t > 0. 

(PI) For every t > 0, the function x i— > v(x,t) is continuously differentiable on Aq. 
(P2) Given e > there exists C > such that \Vv(x,t)\ < C for all (x,t) G 
A t x [0,+ao). 

(PS) Given e > and rj > 0, there exists 5 > such that 

\Vv(y,t) - Vu(x,t)| < J] 

for all x, y G A e satifying \x — y\ < 5 and for all t > 0. 
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The notation Vv(x,t) stands for the derivative of v with respect to x. In order 
to prove this proposition, we will need the following lemma and its corollary. 

Lemma 2.8. Suppose g G L 1 (M d ) is nonnegative, radially symmetric decreasing, 
and supported in B(0, 1). Suppose also that $ G C"(M d \{0}) n L} oc {R d ). Then <f> * g G 
C(K d \{0})nL ; 1 oc (R d ) and 

\\$*9h~(A t )<\\9\\L H B c/ J sup \m\+ ll ®t H * 2) }- ( 2 - 22 ) 

^ e/2<\y\<2 \ a e/2\ > 

Proof. Since $ belongs to Lj oc , by Young's inequality $ * g is also in L] oc . We 
now prove estimate (|2.22p . Fix i in 4 e . 

\($*g)(x)\ < f \${x - y)\ g{y)dy + f \$(x - y)\g{y)dy (2.23) 

< ( sup \\g\\ L i [Bt/2) (2.24) 

+ ( sup g(yj)\\n L HB 2 )- (2.25) 

Since g is radially symmetric decreasing, we have that g(z) > swp yeA g(y) for almost 
every z G B e / 2 - Therefore we obtain 

IMUi(£ e/2 ) > \B e / 2 \ sup g(y) 
y£A c/2 

which, combined with (|2.25|) . leads to the desired estimate. We now prove that <£> * g 
is continuous. Reasoning as above we obtain that if x G A e then 

|(* » g)(x + />)-(*. g)(x)\ < NUmb,,,) 

{ sup |,( 9 UM( S )| + W1M*™). (2.26) 

^ e/2<\y\<2 \- t >e/2\ > 

We conclude that $ * g is continuous using the fact that $ is uniformly continuous 
on compact sets which do not contain the origin and the continuity of the translation 
/l ->■$(• + h) from R d to L 1 (B 2 ). □ 

Recall that P flr ,(R d ) is the space of probability measure fi which can be written 

/i = mS + g 

for some m > and for some nonnegative, radially symmetric decreasing function 
g G L 1 (R d ). From the previous Lemma we directly obtain the following corollary: 

Corollary 2.9. Let K : R d ->■ E 6e a potential such that K XiXj G C(K d \{0}) n 
Lj oc (lR ). TTien £/ie family of functions {K XiXj * /j, : /i 6 PflDfl'')} is uniformly 
bounded and equicontinuous on every annulus A e , e > 0. To be more precise we 
have: 

\\K XiXj * fi\\ L °° (Ae) < sup \K XiXj (y)\+ Kxi *> L1(B2) forall^er RD (R d ) 

e/2<\y\<2 \- t >e/2\ 

(2.27) 
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And also: given e > and rj > 0, there exists S > such that 



- (K XiXj *n){y)\<ri 



(2.28) 



/or all x,y £ A e satisfying \x — y\ < S and for all /i 6 PrdOR ). 

Proof. Since /i G 7-^.0 (R d ) it can be written /i = mS + g where m G [0, 1] and 
5 satisfies the hypothesis of Lemma 12.81 So K XiXj * fi — mK XiXj + K XiXj * g and it 
is easy to conclude using Lemma 12.81 with $ = K XiXj . The second statement is a 
consequence of (|2.26p . □ 

We now prove Proposition 12. 71 

Proof of Proposition \2.7\ Let us first check that v(x,t) is continuous on Aq x 
[0, +oo). Continuity with respect to time comes from the fact that 1 1— > p(t) is narrowly 
continuous together with the fact that the kernel <fi appearing in formula (|2.6[) belongs 
to C°([0, +oo)). Continuity with respect to space comes from the Holder continuity 
of 4>. To prove (Pf), (P2) and (P3), note first that if K{x) = \x\ a \ a > 2 - a, then 
K XiXj G C(K d \{0}) n Lj oc (M. d ), and then use Corollary El □ 

Definition 2.10. V is the space of velocity fields v : -8(0,1) x [0, +oo) — s- R d 
which are radially symmetric and pointing inward (i.e. v(0,t) = and v(x,t) = 
— X(\x\ ,t)x, \x\ > 0. for some nonnegative function A : (0,1) x [0, +oo) — > M.) and 
which satisfies (PO)-(PS). 

Obviously a velocity field defined as in the statement of Proposition ^ . 71 belongs to 
V (the fact that it points inward comes from formula (|2.6[) together with the positivity 
of the kernel (f>). We now investigate properties of flow maps generated by velocity 
fields in V. 

Proposition 2.11. Suppose v G V. Then there exists a unique continuous 
function a : B(0, 1) x [0, +oo) i-> -B(0, 1) satisfying 



Moreover if the point (xo,to) ^ P>(0, 1) x [0, +oo) is such that er(xo,io) ®> then the 
mapping x H- a(x,to) is continuously differentiable at xq and we have 



In particular, if R(t) is such that {a 1 ) 1 ({0}) = -B(0, R(t)), then a 1 is a diffeomor- 
phism from 5(0, l)\B(0,R(t)) to M d \{0}. 

Proof. The global existence and forward uniqueness of solution of the ODE x = 
w(x, t) simply come from the fact that v is Lipschitz continuous with respect to space 
away from the origin (because of (P2)) together with the fact that v is pointing inward. 
Indeed, solutions can be continued as long as they are in B(0, 1)\{0} and since v is 
pointing inward, the only way for a solution to escape this domain is to reach the 
origin. When a solution reaches the origin, it stays there forever in accordance to the 
fact that u(0, t) = for all t > 0. 

The differentiability of cr* on 5(0, 1)\B(0, R(t)) and formula (f230)) are more del- 
icate. In classical ODE textbooks, such results are obtained under the assumption 
that v(x, t) is continuously differentiable in both space and time. In our case v is 
continuously differentiable in space but only continuous in time. However, by revisit- 
ing classical proofs, one can easily check that assumption (P3) is enough to obtained 




for all (x, t) G B(Q, 1) x [0, +oo). 



(2.29) 




(2.30) 
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differentiability of the flow map as well as formula (|2.30[) . This is done in section RTT1 
of the appendix. □ 

Since the flow map a 1 generated by a velocity field in V is a diffcomorphism 
from B(0, 1)\B(0, R(t)) to B(0, 1)\{0}, we can use the change of variable formula to 
express the push forward of a measure in ■pRi)(R d ) by er*: 

Corollary 2.12. Suppose v £ V and let a be the associated flow map provided 
by Proposition 1 2. 1 1\ Let /i = too<5 + go G VRD(R d ) and assume that supp/i G B(0, 1). 
Then 

a t ^L = m{t)8 + g{t) 
where m(t) G M+ and g(t) G L 1 (M d ) satisfy : 

m(t) = too + / 9o(x) dx 



(2.31) 
(2.32) 



2.3. Radially symmetric decreasing profiles are preserved. 

In this subsection we give a heuristic argument (which will be made rigorous 
in the next subsection) explaining why radially symmetric decreasing profiles are 
preserved by the aggregation equation when 2 — d < a < 2. We recall here that it is 
observed numerically that when a > 2. radially symmetric decreasing profiles are not 
preserved. A key ingredient in our argument is the known fact that the convolution 
of two radially symmetric decreasing functions is still radially symmetric decreasing 
(see [54] for example). For completeness we give a quick proof of this fact: 

Lemma 2.13. Suppose g G L 1 (R d ) with compact support and f <G Lj oc (R d ). If 
f and g are nonnegative radially symmetric decreasing functions, then f * g is also a 
nonnegative radially symmetric decreasing function. 

Proof. The mononicity of g allows us to use a "layer cake" decomposition of g, 
namely 

POO 

g(x) = / XB(o,f{g)){x)dg 
Jo 

where r(g) denotes the inverse function of g(r) and Xb(o,s) denotes the characteristic 
function of the ball of radius s. Thus 

/•OO 

f*g(x)= f * XB(o,f( g )){x)dg (2.33) 
Jo 

and we note that the integrand of (|2.33|) is monotone decreasing because the char- 
acteristic function of a ball convolved with a nonnegative Lj oc monotone decreasing 
function is itself monotone decreasing. By integrating a monotone integrand with 
respect to g we obtain the monotonicity result for / * g. □ 

We now present the heuristic argument. Let us assume that u(x, t) is a smooth 
solution of the aggregation equation. Clearly we have: 

u 

— + Vit • v = (AK * u)u. (2.34) 
ot 
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Suppose that AK is locally integrablc, nonncgative and radially symmetric decreasing. 
When VK(x) — x\x\ a ~ 2 , these hold true if and only if 2 — d < a < 2. We then use 
Lemma 12.131 to see that if for some to > 0, u(-,to) is radially symmetric decreasing 
then the right hand side of (|2.34[) is also radially symmetric decreasing at to- This 
indicates that the rate of change along the characteristic is greater the closer we 
are to the origin. Therefore the solution is expected to remain radially symmetric 
decreasing for t > to. For the special case of the Newtonian potential, AK * u = u 
and monotonicity is similarly preserved - this is discussed in more detail in Section 3J 

2.4. Proof of Theorem [HB 

The proof is inspired by the work in [28| . where global existence of measure 
solutions for some kinetics model was obtained by using a fixed point iteration in the 
space of probability measures endowed with the Wasserstein distance. 

Let VK(x) = x\x\ a ~ 2 , a G (2 - d, 2), and let p init G V RD (R d ) with supp(p lmt ) C 
5(0,1). Define: 

Pa(t) = Pinit Vt £ [0, +oo) 

j-(VK* Po (t))(x), ifx^O w , crn ^ , 
vo[x,t) = i Vte[0,+oo) 
10, if x = 

erg : K d — > K d := flow map associated with vq 
and for n > 1 define recursively 

pn{t) = <J t n _ l #p lmt Vt G [0, +oo) 

v n {x,t) = i Vte[0,+oo) 
10, if x = 

a f n : W d —> M. d := flow map associated with v n . 
Proposition 2.14. 

(i) For all n > 0, p n G C([0, +oo), V aD Q& d )) and supp(p n (t)) C B(0, 1) for all 
t > 0. 

(ii) Given e > 0, there exists L e > such that 

\v n (x,t) - v n (y,t)\ <L e \x- y\ 

for all x, y G A t , for all t > 0, and for all n > 0. 
(Hi) There exists a constant 9 G (0, 1] depending only on a such that the following 
holds: Given e > ; there exists C' e > and 8 > such that \t — s\ < 5 implies 

\v n (x,t) - v n (x,s)\ < C E \s - tf 

for all x G A e and for all n > 0. 
(iv) Pn+i (t) y Pn(t) for all n > and all t G [0, +oo). This implies \v n +i{x, t)| > 
\v n (x,t)\ for all (x,t) G R d x [0, +oo) and for all n>0. 

Before we prove this proposition let us explain how it will be used in the proof of 
Theorem 12.11 Because of statements (ii), (iii), (iv) and the bound |w„(a;,t)| < |a;| Q_1 
(see (|2.9p ). we can use the Arzela-Ascoli theorem to conclude that the v^s converge 
uniformly on A t x [0, +oo) to some function v which is Lipschitz continuous in space 
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and Holder continuous in time, with same constants L e and C e . Since e can be 
chosen as small as we want, v(x,t) is well define on 5(0, 1)\{0} x [0, +00). Let 
v(0,t) — so that v is now well defined on 5(0, 1) x [0, +00). This velocity field 
v(x, t) generates a flow map <r* : 5(0, 1) — > B(0, 1) and from this flow map we can 
construct p{t) = o- t #Pi n it- In Proposition 12.151 it will be shown that er„ converges 
uniformly to a on 5(0, 1) x [0, 00). This implies in particular that for a given t, p n (t) 
converges narrowly to p(t). The narrow convergence preserves the monotonicity (see 
Proposition [HUl of the Appendix), and therefore p(t) is radially symmetric decreasing. 
In order to prove that the radially symmetric decreasing function pit) and the flow 
map cr* obtained by the above limiting process satisfy (|2.1[) and (|2.2[) we just need to 
show that v(x,t) = — (Vif * p(t))(x) for x ^ 0, and this fact will follow easily from 
passing to the limit in the relationship v n (x,t) = —{S/K * p n (t))(x) for x ^ 0. 

Proof of Proposition \2.14\ Let us first prove (i) . The initial iterate po it) = Pinit 
obviously belongs to C([0, +00), Viiz>(M d )) with supp(p (i)) C 5(0,1) for all t > 0. 
Assume that p n g C([0, +00), V RD (R d )) with supp(p„(i)) C 5(0,1) for all * > 0. 
From Proposition 12 . 71 v„ E V and from Corollary 12.121 

Pn+i(t) = m n+1 (t)6 + g n+ i (t) (2.35) 
m n+1 (t) = to + / g (x) dx (2.36) 

Sn+i(*,t) = (^T^T i^r 1 ) W fOT x * 0. (2.37) 

Here to and go are such that /9j„it = m S + g - Also from Proposition 12 . 1 ll we know 
that detVcr^ satisfies 

detVcr^(x) = cxp f (div v n )(a^(x), s)ds (2.38) 
Jo 

for all (x,t) such that cr^x) ^ 0. Since we have assumed that p n (t) is in VnD(^. d ) 
with compact support, and since for a G (2 — d, 2) Aif is nonnegative, radially 
symmetric decreasing and locally integrable, we know from Lemma 12.131 that the 
function x 1— > — divu„(a;,t) = [Ai^ * p n (i)](x) is nonnegative, radially symmetric and 
decreasing. Since |a;| < \y\ implies |er*(x)| < |cr^(y)| one can easily see from (|2.38p 
that 

x H» ttt i s nonnegative, radially symmetric and decreasing. 

det Vct^(x) 

Then we easily see from f|2 . 3T[) that, since go is radially symmetric and decreasing, so 
is x i-» g n+ i(x,t). 

Let us now remark that the estimate |u n (:e, i)| < \x\ a 1 together with Lemma l6.5l 
of the Appendix lead to the following: if a 6 (2 — d, 1) then 

\ai(x) - a s n (x)\ <C a \t-s\^ (2.39) 

for all x € B(0, 1) and for all t, s > 0. Here C a := (2 - a)^. If a £ [1,2) then 
v {x, t) < 1 on 5(0, 1) x [0, +00) and therefore we get 

|<T*(i)-<(x)|<|t-a| (2.40) 
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for all x G -8(0, 1) and for all t, s > 0. Using Lemma 15751 from the Appendix, together 
with (|2.39l) we obtain that, if a G (2 — d, 1) then 



W 2 {p n+1 {t), Pn+1 {s)) < ||<7* -O-^IUo^o,!)) < C Q |t - s|*-«. 

This prove that i i— >• p„ +1 (t) is Holder continuous with respect to when a G 
(2 — d, 1). If a G [1, 2), we obtain from (|2.40p that t n- p„ + i(t) is Lipschitz continuous 
with respect to W2. 

Statement (ii) is a direct consequence of Corollary |2.9l We now prove (iii). Sup- 
pose 2 — d < a < 1. Recall from Lemma T2.6I that <fi is 7-Holdcr continuous for some 

1 

7 G (0, 1]. Choose 8 such that \t — s\ < 8 implies C a \t — s\ 2 - a < e/2. Using Lemma 
and estimate (|2.39[) we obtain that \t — s\ < 8 and x £ A e implies that 



\v n (x,t) - v n (x, s)\ = \x\ a 1 







. dpinit(r) (2.41) 



r+00 

^^r 1 - 1 cla^^-a^MUpinuir) (2.42) 





<cC Q |xr- 1 - 7 |i-s|^. (2.43) 

The first equality is a simple consequence of formula (|2.6[) . the fact that p n (t) = 
^n-i^Pinit and the definition of the push forward. Note that C a \t — sl 5 ^ < e/2 
and (|239)1 imply that ^ui^ - f %p^ < V 2 for £ e A e . This allowed us to use 
Lemma [2.61 in order to go from (|2.41[) to (|2.42[) . The case a G [1,2) is dealt with 
similarly. 

We finally prove (iv). Obviously pi{t) >- po{t) = pinit for all t > 0. Assume that 
for a given n, p„(t) >- p n -i{t) for all < > 0. Then (|2.8p implies |u n (x,t)| > |u„_i(x, t)|. 
Lemma 13.21 which is proven in the next section, implies then that p n +\{t) >~ p n {t) for 
all t > 0. □ 

As already mentioned, (ii) (iii) and (iv) imply that the sequence {v n } converges 
uniformly on A e x [0, +00) to some function v (which is Lipschitz continuous in space 
away from the origin). Setting v(0,t) =0wc obtain a velocity field well defined on 
5(0,1) x [0, +00). This velocity field v(x,t) generates a flow map a* : 5(0,1) — > 
5(0,1). 

Proposition 2.15. a n (x,t) converges uniformly to a(x,t) on 5(0, 1) x [0, +00). 

Proof. Let e > be fixed. From formula (|2.6j) it is clear that |i>o| is strictly positive 
away from the origin. Since l^n+il > \v n \ we have that \v\ is also strictly positive away 
from the origin. Therefore there exists a time T e > such that a Te (5(0, 1)) C 5(0, e). 
Choose N so that n> N implies \\v — v n\\ ^ ^A^x[o.T e ]) 

Case 1: Assume first that (x, t) G 5(0, 1) x [0, +00) is such that |cr* (a?) | > e. Note 

15 



that such a t is necessarily smaller than T e . For all r < t and for all n > we have 



|a T (x)-<(x)| < / «)| ds 

Jo 

< / |i>(cr s (;r),s) - v n (a s (x),s)\ + \v n (cr s (x), s) - v n {cr s n (x) , s)\ ds 
Jo 

< T \\ V ~ V n\\ L ^ {Ae x[^T])+ L ^ / W 8 (x)-CT^(x)\ds. 

J 

We have use the fact that |<r*(x)| > e implies that |c s (x)| > e for all s < r < t. We 
have also use the fact that, since \v\ > \v n \, cr^(x)| > \a s (x)\ > e for all s < r < t 
and for all n > 0. Using Gronwall's lemma and the fact that t < T e we obtain that 
for n > N: 

\a\x) - al(x)\ < T e \\v - «„|[ £ . (x . x[0jT , ]) e L = T = < e. (2.44) 

Case 2: Assume that (x,t) G -6(0, 1) x [0, +oo) is such that \x\ < e. Since the 
velocity fields v and u„ are focussing we clearly have that |cr*(x) — cr^(x)| < 2e for all 
n. 

Case 3: Assume finally that (x,t) £ B(0,1) x [0,+oo) is such that |cr*(x)| < e 
and |x| > e. Since r i— > cr T (x) is continuous there exists a time s £ [0,t] such that 
|cr s (x)| = e. So from case 1 we get |cr s (x) — <r^(x)\ < e for n > N. Since |cr s (x)| = e 
we have that |er*(x)| < 2e for n> N. Since s<twe have |cr^(x)| < |c*(a;)| < 2e for 
n > N. Therefore |cr*(x) - cr^(x)| < 3e for all n > N. 



We are now ready to prove the Theorem 12. II 

Proof of Theorem \2.1\ Define p(t) := o- t #Pinit- Recall that from Lemma T6.3I of 
the Appendix 

W 2 (p,Pn)~ sup W 2 (p(t),p n (t)) <\\a-a n \\ LX{Bix[0+oo)) . 

tG[0,oo) 

So from Proposition 12 . 1 51 we get that W2(pn — p) — > 0. This implies in particular that 
for every t £ [0, +oo), p n {t) converges narrowly to p(t). Since narrow convergence 
preserves the monotonicity fLcmma l6.4l of the Appendix), we know that p(t) is radially 
symmetric decreasing. 

We are now going to prove that p and a satisfy (|2.1[) and f|2.2f) . Since p is defined 
by p(t) = o- t ^pi n it where a* : M. d — > K d is the flow map associated to the velocity field 
v(x, t), we just need to prove that v(x, t) = —{VK * p(t))(x) for x ^ 0. This is obtain 
by passing to the limit in the relation v n (x,t) = —(VK * p n (t))(x) for i^0. Indeed 
v n converges pointwise to v in A$ x [0,+oo). And since for fixed t, p n (t) converges 
narrowly to p(t), we obtain from (|2.6p that VK * p n converges pointwise to VK * p 
in A a x [0, +oo). □ 

3. Instantaneous mass concentration. 

The recent work of [12j [33] concerns local well-posedness of the problem with 
initial data in L p . One can prove a sharp condition on p for local well-posedness by 
considering a family of initial data that behave as a powerlaw near the origin. Such 
initial conditions satisfy the monotonicity assumptions considered in this paper. In 
this section, using existence results from the prior section and a bootstrap argument, 
we prove results about the behavior of these solutions as measure solutions that 
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concentrate mass. Such results are not discussed in the prior literature for the singular 
power law potential K (x) = \x\ a , a < 1. 

More specifically in |12j it was proven that the aggregation equation with potential 
VK(x) = x \x\ a ~ , 2 — d < a < 2, is locally well posed in any L p -space with p > 



Note that given (i <E 2 . 1) the function 

h{x) 




d+a-2 " 

if \x\ < 1 
otherwise 

belongs to the critical space L d +<*- 2 (R d ) but does not belong to any LP space with 
p > d+ ^_ 2 . In |33| it was proved that if the initial data is exactly equal to h(x) then a 

solution of the aggregation equation instantaneously leaves the space L d +°- 2 . In this 

d 

section we go a little further and show that the solution not only leaves L d + a ~ 2 but 
also instantaneously concentrates some point mass at the origin. Our results make use 
of the existence theory from the previous section. Also compared to the work in |12) 
and |33j . our argument here is local in essence and holds for any radially symmetric 
decreasing initial data which is locally more singular than h(x) at the origin. The 
main theorem of the section is the following: 

Theorem 3.1. Let VK(x) = x \x\ a ~ 2 , 2 - d < a < 2. Suppose p mlt G P RD (R d ) 
is compactly supported and absolutely continuous with respect to the Lebesgue measure. 
Suppose that there exists c > 0, r$ > and fi G ( rf+ rf~ 2 1 1) such that the density Ui n n 
of Pinit satisfies 

c 1 

Suppose finally that p 6 C([0, +oo), T > RD0& d )) satisfies the Lagrangian formulation 
(|2.1[) - (|2.2p of the aggregation equation. Then p(t)({0}) > for all t > 0. 

3.1. Comparison principles. 

In this subsection we derive a few comparison principles which will be necessary 
in order to make the arguments local. 

Lemma 3.2. Suppose Ui,i>2 € V and \v\ \ > \v2\- Then 

<t\#H >- o-2#M for all fj, e ^(K d ) and t > 

where o\ and 02 are the flow maps associated to V\ and vi respectively. 

Proof. Since v-i G V the flow map a\ is invertible away from the origin. Define 
A|(a;) = {o~2)~ 1 {x) if x ^ and A|(0) = 0. One can then easily check that 

(a t 1 oA*)#(a*# Ai ) = ^# M 

Moreover since \v±\ > \v2\ we have that |(crJoA|)(a;) < |x|, which concludes the 
proof. □ 

Lemma 3.3. Suppose »£V. Suppose also that p, v € V]i(K d ) and p>~ v. Then 

a f ifp y a l jfv for all t > 

where a is the flow maps associated to v. 

Proof. Since p >- v there is a map P satisfying |P(x)| < |x| such that p = P^v. 
As in the previous lemma, define A*(x) = (a*)~ 1 (x) if x 7^ and A*(0) = 0. One can 
then easily check that 



(a t oPoA t )#(a t #p)=a t ^ 
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and \(a f o P o A t )(x)\ < \x\ which conclude the proof. □ 

The following definition will be needed in order to compare two measures of 
different mass. 

Definition 3.4. Suppose p e V R {M. d ) and p, e M R {R d ), with p(R d ) < 1. We 
write p\> p if there exists a measure v G V R (M. d ) such that 

p>v and v{A) > p(A) MA G B(R d ). 

In view of (|2.7[) and (|2 .8[) it is clear that: 

p>p => \VK * p\> \VK * p\ (3.2) 

The following Lemma will be useful in order to make localized comparisons. 

Lemma 3.5. Suppose v\,V2 G V and \v\\ > 1 1^2 1 in B(0,2R) x [0, +00) . Suppose 
also that p G V R (R d ), p G M R (R d ) and p>p. Then 

<A#p > o-2#(hXb(om)) forallt>0, 

where o\ and o 2 are the flow maps associated to V\ and V2 respectively, and Xb(o,R) 
is the indicator function of the set B(0,R). 

Proof. Since p > p there exists a probability measure v such that p >- v > fx. 
Let £(x) be a smooth radially symmetric function which satisfies £(x) = 1 if x < -R, 
= if \x\ > 2R and x( x ) < 1 for all x G K d . The velocity field v 3 (x,t) := 
W2(x,i)£(x) is still in V. Moreover we have ^3 1 < \v±\ for all x G M. d and t > 0. We 
can therefore use the two previous Lemmas to obtain that 

<A #P >- o-l#P y crl#v > a\#p > <t\#{}iXb%r\) 

The last two inequalities are a simple consequence of the definition of the push- forward 
together with the fact that v > p > PXb[o,R\- Finally, note that since U3 = v-i on 
B(0,R) x [0,+co), then 4#G« Xb[o,r]) = ^2#{pXb^ : r])- □ 

3.2. Proof of Theorem 13.11 by bootstrap argument. 

Fix a G (2 — d, 2) and define the functions 

fe,ro( x ) = , ,d+I-2+e X-B(0,r o ) 0*0 for 6 G (0,1) (3.3) 

\x\ 

9r ( X )= , t d+ a -2 XB(0,r )(x) (3.4) 

\x\ 

h P,M = (-\ n \x\)P XB M^ fOT 13 e ( 2 ' ^' (3 ' 5) 

Note that at the origin f e , ra is more singular than g ro which itself is more singu- 
lar than h/3 t r . In |12j it was proved that if a = 1 and the initial data is exactly 
equal to Cf t ^ ro (x) (C is a normalizing constant) then a Dirac delta function appears 
instantaneously in the solution. The proof relied on the fact that solutions of the 
ODE x = — (VK * / e)ro )(x) reach the origin in finite time. However this strategy 
does not work with g ro and / e , ro , because solutions of x = — (VK * g ro )(x) and 
x = —(VK * h/3 tro )(x) do not reach the origin in finite time. For that reason we 
will use a bootstrap argument to prove that a delta function appears instantaneously 
when the initial data is equal to or more singular than /is, ro G L d + a ~ 2 (R d ). Roughly 
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speaking, wc will show that the velocity field — V-ftT * hp ro instantaneously deforms 
ha,ro into a function more singular than g ro , then we will show that the velocity field 
— VA" * g ro instantaneously deforms g ro into a function more singular than f e , ro , and 
finally wc will use the argument from [12] to show that the velocity field —VK * f eyro 
deforms f tt r in such a way that a delta function appears instantly. 
The following definition is consistent with Definition 12.21 

Definition 3.6. Given a radially symmetric, non-negative function u G L 1 (W i ), 
we define u G L ((0, +oo)) to be the unique function satisfying 

/ u(r)dr = / u{x)dx for all r*i,r2 > 0. 

J n Jr 1 <\x\<r 2 

In other words, u(r) = w(r)cjd? ,<i_1 . 
With this notation we have: 

fe,r {r)=U> d ra \ +e X[0,r o )( r ) ( 3 - 6 ) 
9r (r) = UJ d ^ZlX[0,r o ]W (3.7) 
h,r (r)=Ud r «-i (_ lnT .)/3 X[o,ro](r). (3.8) 

We remind the reader that by Lagrangian solution we mean a function p(x, t) that 
satisfies the Lagrangian formulation (|2.1j) - (|2.2[) of the aggregation equation. 

Proposition 3.7. Let p G C([0, +oo) 1 VB,D(M d )) be a Lagrangian solution of 
the aggregation equation with compactly supported initial data pmit and potential K 
satisfying VK(x) = x \x\ a ~ 2 , 2 — d < a < 2. Let v(x,t) = —(VK * p(t))(x) be the 
associated velocity field. 

(i) If pinit > c/ e ,r for some c, ro > and e G (0, 1), then there exist R,C > 
such that 

\v(x,t)\ > C |x| 1_e for all (x,t) G 5(0, R) x [0 + oo). 
(ii) If pmit » c g ro for some c, r > 0, then there exist R,C > such that 
\v{x,t)\ > C \x\ (— In for all (x,t) G B(0, R) x [0 + oo). 

(Hi) If pinit>ch/3.r for some c, ro > and j3 G ( ~*~^~ 2 , 1), then there exist R,C > 

\v{x,t)\ > C \x\ {-\a\x\) l - f3 for all (x,t) G B(0,R) x [0 + oo). 

Proof. Let us prove (i). On one hand, from fj3 . 2|) we see that \v(x, 0)| > 
c \(VK * f e . ro )(x)\ for all x G M d . On the other hand, since the velocity field is 
always pointing inward (this is due to the positivity of (/)), we have that p{t) y p(0) 
for all t > 0, and therefore from (|2.8p we get that > |u(x, 0)| for all x G K d 

and t > 0. So we only need to show that |(Vif * / c .r )( :z; )l ^ C M l e in some neigh- 
borhood of the origin, and this estimate follows easily from Lemma 12.41 Indeed, by 
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Lemma 12.41 wc have for |x| < r 



. a-1 

r \ l \x\ 



■|x| 

\VK * f e , ro (x)\ =Lo d I ^(nlP) r-'dr (3.9) 



x 7 \ r 



o i r \ I r 



^'U/IhAhJ * '* (3 ' 10) 

>w d Ci / r~ £ dr + cj d C 2 |:r| / r^^dr (3.11) 
Jo J\x\ 

>uj d C^ (3.12) 

1 — e 

where C\ = fnf[ 0i i] 4> — 0(1) and C 2 = inf( 1+00 ) <j){r)r 2 ~ a > 0. 

Let us now prove (ii). Reasoning as above we see that it is enough to show 

that |(VA' * g ro )(x)\ > C \x\ (— ln|x|) in some neighborhood of the origin. Then the 
argument is similar. From p. lip with e = we get 



\VK*g ro {x)\ >WdC 2 |x|ln 



which yields to the desired estimate. 

To prove (iii) it is enough to show \(VK * hp iro )(x)\ > C \x\ (— In \x\Y~ 13 in some 
neighborhood of the origin, and the argument is similar. In this case we have 

f r ° 1 dr 

\VK * uq(x)\ > u>dC2\x\ 



I logr|. 



which yields to the desired estimate. This last estimate was derived independently in 
[33]. □ 

The ODE's 

r = -Cr 1 -*, f = -Cr{- In r) and r = -Cr{- In r) 1_/3 

suggested by the previous proposition have explicit solutions and their flow maps are 
respectively: 

t,, / \ i (r e - eCt) 1 /' if r > (eCt) 1/e , n . 

aUr) =ai(r,t) = I v ; ) (3.13) 

1W V ' ' [0 if r < (eCt) 1 /' V ' 

4(r)=V2(r,t) = r e0t (3.14) 

4(r)=* 3 (r,t) = e-( C ^^r) 1/ * (3.15) 

Solutions of the first ODE reach the origin in finite time but solutions of the other 
two ODE's only approach the origin as t -> oo. Corresponding to the flow maps 
tii : [0, +oo) x [0, +oo) -)• [0, +oo) there arc flow maps 5, : R d x [0, +oo) R d 
defined by Si(x,t) = crMxl , t)-^. The Si are the flow maps associated to the 

velocity fields Wi(x) = —C\x\ ~ e rfr, W2(x) = — C \x\ (— In xQjfy, and 103(0:) = 
— C (— In |x|) 1 ~ ,3 |^| . Let u £ L 1 (R d ) be a radially symmetric, non-negative func- 
tion. It is clear from (|3.13[) that S{#u has a point mass at the origin if u has non-zero 
mass in B(0, (eCi) 1 / 6 ). On the other hand, because St^ and S\ are smooth invertiblc 
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maps, Syftu and S^^u are continuous with respect to the Lebesgue measure, and by 
the change of variable formula, we have 



(Sl#uT(r) = {oUu){r) = u(r/(r)) i = 2,3 (3.16) 

where r*(r) = (a*)" 1 ^) (3.17) 

PROPOSITION 3.8 (Bootstrap). Let p G C([0, +oo), V RD (R d )) be a Lagrangian 
solution of the aggregation equation with compactly supported initial data pinit cind 
potential K satisfying VK(x) = x\x\ a 2 , 2 — d < a < 2. 

(i) If pi n it > c/ e ro for some c,vq > and e G (0, 1), then p(t)({0}) > for all 
t > 0. 

(ii) If pinit > c g ro for some c,ro > 0, then for any t > there exist constants 

Ci, ri > and e G (0, 1), such that p(t) > c\ f e . ri . 
(Hi) If pinit > c hp yro , for some c, ro > and (3 G ( rf+ ^~ 2 , 1), i/ien /or any £ > 

there exists constants ci,r\ > such that p{t) t> c\ g ri . 
Proof. Let us prove (i). Let S[(x) = Si(x,t) be the flow map generated by the 
velocity field w\(x) = — C\x\ _e A suggested by Proposition 13.71 From Lemma [3?5l 
and Proposition 13.71 we then obtain that 

p(t)>5*#c/ e>ri (3.18) 

for r\ small enough and for all t > 0. Let us fix a t > 0. Since / e ri has non-zero 
mass in B(0, (eCt) 1 ^), it is clear from (|3.13[) that the measure S\#cf t ^ ri has a point 
mass at the origin. Then by (|3 . 1 8[) we conclude that p{t) also has a point mass at the 
origin. 

Let us now prove (ii). Again Lemma [3.51 and Proposition 13.71 imply that p(t) > 
S^i^cOn for r\ small enough. As already mentioned S\^cg ri is continuous with 
respect to Lebesgue measure. We are going to show that given any t > 0, S\j^g ri > 
C2fr 2 ,e for some constant C2, r^ > and e G (0, 1) which will conclude the proof of (ii). 
Let T-|(r) = {o~2)~ l { r ) = ?" e where o- 2 (r) is defined by (|3.14|) . Using the change of 
variable formula, we get that 

^2*9r x ){r)=g r A4{r)) ^(r) 

Ud -ct„e ' 



( r ) 



e r e for r small enough 



,,a-l + (2-Q)(l-e- ct ) 



Since 2 — a > it is clear that (o"2#<?r-i)(?') > C2.f r2 ^(r) for r small enough. 

Let us now prove (hi). Once more Lemma 13.51 and Proposition 13.71 imply that 
p(t) f> S l 3#cft. ( g. ri for r\ small enough. Let us fix t > and show that Sl#chp iri > C2<? r2 
for ri small enough. In view of (|3.7j) it is enough to prove that 



limr"- 1 (^#V ri (r)) > 0. (3.19) 
Let t\(t) = (<J^)~ 1 (r) and note that 

/ / IN A 17 " 

(3.20) 




From now on we drop the lower subscript. From the change of variable formula we 
have 

a t #h(r) = h(r t (r))^ r (r) (3.21) 

t( U-a Q t 

-(r) (3.22) 



i \P dr 



where we have used ((3~2T)|) to go from (j3~2"2"j) to (j3~2"3")) . Then note that using (j3~2T)]) 
again we get 



T*(r) dr \T t (r) J \ \ r 

which combined with p.23|) gives 



r / r 



-. ( ^ W ).(iM)-^ i+ (,„i)y- 2 ( ln i)- (3 , 4) 

t (r) \ / cpi \ / 1 1 



i-TT^I ( ln r) ( 3 - 25 ) 



r 



1 /r*(r)\ 2 ~ Q / 
> — In - for r small enough (3.26) 



2 \ r J \ r 
Using (|3.20j) and doing a Taylor expansion we find that 

i / iv-* 3 

> —ct I In — J for r small enough (3.29) 



Combining (j3~26| and (pT2"9j) we get 



| > ln(l/2) + ^ (2 -a)ct fin-") ^-pTnln- 



for r small enough. Since 2— a > it is clear that linv^oo In yr a 1 ya t ^h(r)J J = +oo 

which implies (|3.19[) . □ 

We now prove Theorem 13. II 

Proof of Theorem \3.1l If pmit > c ft^, ro f° r some for some c, ro > and /3 € 
( j~ 2 , 1), we can apply the previous proposition to get that for any t\ > 0, p{t\) t> 
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c g ro for some different constants c, r > 0. Applying the proposition again, we 
get that for any i 2 > ti , p(p2) > c /em f° r some other constants c, ?* > and for 
some e € (0, 1). Applying the proposition one last time we get that for any > 
P(*3)({0}) > 0- Since ti < t 2 < t 3 can be chosen arbitrarily small, this conclude the 
proof. □ 

4. Newtonian potential case. 

An even more singular case is that of the Newtonian potential, a = 2 — d in 
general, with K(x) = log \x\ in the special case of 2D. Without loss of generality we 
use the normalization for K that yields A (if * p) — p, i.e. the fundamental solution 
of the Poisson equation. This simple fact localizes the dynamics as compared to 
the nonlocal case studied in previous sections. In Eulerian coordinates, for smooth 
densities, we have 

p t + v ■ Vp = p 2 . (4.1) 
Recall that for radially symmetric problems, the Laplace operator is 

Likewise we have the following formulae for the gradient and divergence operators: 

or 

where r is the unit outward pointing radial vector and 

A' 1 9 d-1 

dlV V = — j—r — — V V. 

r a 1 or 

Using the latter formula we can rewrite v above in terms of p simply by inverting 

divv = —p: 

<r) = -^ x [s^p{s)ds:=-^ (4.2) 

where m(r) is proportional to the mass contained inside a ball of radius r. Thus it 
makes sense to rewrite the evolution equation (|4.1[) in mass coordinates - in general 
regardless of the kernel it is 

m t + vm r = 0. (4-3) 
However this greatly simplifies in the special Newtonian case. Formula (|4.2p gives 

m * - — = °" 

By changing variables to z coordinates, where z = we have the inviscid Burgers 
equation, 

m t — mm z = 0. (4-4) 

The transformation to equation (|4.3[) is well-known; the transformation to the z 
variable appeared in |14j in the context of a viscous version of our problem arising 
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in astrophysics. Here we use the classical conservation law theory for the inviscid 
(purely transport) problem to prove that monotonicity is preserved by the flow in 
all dimensions. In one dimension it is known that K = \x\ can be transformed 
to the inviscid Burgers problem see e.g. |18) . The connection to Burgers equation 
allows us to prove quite a lot about radially symmetric solutions of the aggregation 
equation with Newtonian potential, by directly connecting to the classical theory of 
conservation laws. We consider three cases: (a) monotone decreasing radial densities 
for which we have a unique forward time solution; (b) general radial densities for 
which we have existence of solutions but uniqueness requires the specification of a 
jump condition (akin to choosing a particular entropy-flux pair for the definition of 
distribution solution, and; (c) the case of radially symmetric signed measures for 
which one requires an additional entropy condition in order to have a unique solution. 
All of these cases can be distinguished by the known properties of the inviscid Burgers 
equation |4] |4l?l |57] . 

4.1. Case 1: p € Vrd(^- ) - existence of unique classical solutions. 

In the case of radially symmetric monotone decreasing probability measures, we 
have unique classical solutions by virtue of the fact that the corresponding flow field 
v is Lipschitz for r > 0. Monotonicity is preserved by virtue of the localization of the 
equations as described above. More specifically, we have the heuristic that p satisfies 
pt = p 2 along characteristics so the initial ordering of the density is preserved provided 
that the characteristics remain well ordered and are well defined. We can prove this 
to be the case by going to the mass coordinate formulation above. The condition 
that the characteristics remain well defined is akin to proving that shocks will not 
form from any initial data satisfying the monotonicity condition. If a shock forms - 
which wc define as a singularity in m z in the mass equation (|4.4[) , the first time of 
formation will occur at t s hock = 1/ su P z { m 'init( z )}- So we need the characteristic to 
reach the origin before this time occurs. Denote by z s the location at time zero of 
this characteristic. Then our condition on the shock occurring after the characteristic 
crosses zero is 



since both rrn n u and m' init are nonnegative. Using the definition of the mass m and 
converting back to regular radial coordinates, the above is equivalent to the following 
condition on the density p: 



for all R, which is true for monotone decreasing initial data pinit- The special case of 
the equals sign in (j4.5[) corresponds to the shock happening right when the character- 
istic reaches the origin. There are exact solutions that satisfy this - corresponding to 
a density that is the characteristic function of a collapsing ball. The corresponding so- 
lution in (m, z) coordinates is the well-known Burgers solution of the form —zj (I — t) 
that forms a shock in finite time in which all of the characteristics on an interval 
collapse at the origin simultaneously. This example is the most singular case of the 
general class of solutions considered in this subsection. Since the only shocks that form 
occur at the origin, which is a boundary of the domain, this results in a global-in-time 
classical solution of (|4.4[) for any initial condition mi n it(z) arising from a probability 
density pmAt £ Vrd (K ). The classical solution of the inviscid Burgers equation easily 




z. 



1 




(4.5) 
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gives us a unique solution of the Lagrangian formulation of the problem as well. We 
state these results below: 

Theorem 4.1. Given compactly supported initial data Pi n a 6 'Pfl£>(R d ), define 
minit = Jo s d ~ 1 dp. Then there exists a unique classical solution to equation fi4-4j) 
on the half space (x,t) € (0, oo) x [0,oo) and a corresponding unique solution of 
the Lagrangian mapping formulation of the density transport problem. The solution 
retains its monotonicity property for all time. 

As we show below the situation is much more complicated for general radially 
solutions without the monotonicity condition. Some observations can be immediately 
made using classical results from conservation laws. Moreover these results connect 
directly to related problems in fluid dynamics such as vortex sheet solutions of the 2D 
Euler equations. The next two subsections provide a discussion of these observations. 

4.2. Case 2: p 6 Vfi(M. d ) - existence of unique solutions with jump con- 
dition. 

In the case of general radially symmetric probability densities, we no longer have 
classical solutions. Let us consider the simplest example of data that violates the 
monotonicity condition - that of a uniform delta-concentration on the boundary of 
the ball of radius i?*. Following the mass coordinates, we see that this example has a 
jump discontinuity in m(z) at z* = (R*) d /d. Since to is the characteristic speed, this 
results in a jump in the velocity across the delta-ring. One way to define the solution 
is to consider a distribution solution of (|4.4[) in which case the speed of the shock 
(velocity of the delta-ring) is defined, in z coordinates as s 2 _ s ^ oc j. = (mi + m2)/2, 
i.e. the Rankine-Hugoniot condition associated with equation (j4.4[) . As is well-known 
for scalar conservation laws, we could transform equation (|4.4p by multiplying by any 
function of to, 

(F(m))t - {G{m)) z = 0, F'(m) = f(m), G'{m) =mf(m) (4.6) 

for some function /, yieldling a different jump condition in the weak-distribution form 
of gU), 

_ [F(m)} 

$z — shock — xs-if \ i 7 

[G(m)\ 

where [ ] denotes the jump across the shock. 

By virtue of well-known results for scalar conservation laws, we obtain families of 
weak solutions for the general radially symmetric problem. For a given formulation 
of the form (|4.6|) there exists a unique distribution solution. Uniqueness for inviscid 
Burgers often requires an additional entropy condition. In the case of formulation 
(|4.4p . the entropy condition is automatically satisfied by the monotonicity of to, which 
is guaranteed for any radial probability density, not necessarily monotone. The full 
entropy condition would only be required in the case of non-monotone m such as 
would arise in the case of a signed measure p. 

It would be interesting to know whether there is an optimal choice of shock speeds 
for these under-determined problems. For example one might also consider an optimal 
transport framework in which the best choice of shock speed would be one in which the 
interaction energy is most quickly dissipated. For the aggregation problem this would 
result in the fastest speed possible for the delta-ring which would satisfy a Lagrangian 
formulation of the problem but perhaps not a classical distribution solution in Eulcrian 
variables - even in the m — z framework described above. We note that the entropy 
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solution discussed above, in which the speed of the shock is chosen to be the average 
of the speeds on either side, is a natural generalization of the choice conventionally 
made for 2D vortex sheets, in which v = V K * p rather than v = VK * p, and p is the 
vorticity. For that problem one ascribes a velocity to the sheet that is the arithmetic 
average of the speeds on either side |57j . The frozen time calculation can be made by 
analogy to the incompressible flow problem, however the ensuing dynamics is quite 
different. For the vortex sheet problem the flow is tangential to the sheet so the issue 
of shocks does not arise. For the aggregation problem the flow is normal to the sheet 
and affects the solution on either side of it, because the speed of the shock determines 
the rate at which characteristics on either side of the discontinuity are absorbed and 
the rate at which information is lost in the discontinuity. To summarize, if we define 
a solution as satisfying an equation of the form (|4.6I) in the sense of distributions, 
then we expect a unique solution, however, in the case of jump discontinuities in m, 
the shock speed will depend on the choice of entropy-flux pair as discussed above. 
Moreover we believe even more general examples may exist that could satisfy an 
optimality condition associated with dissipation of the interaction energy. We finally 
briefly mention the case of signed measures below. 

4.3. Case 3: signed measures. 

The case of signed measures introduces yet another source of nonuniqucness of 
solutions, which we briefly discuss. A signed measure corresponds to a non-montonc 
(but solution of the inviscid Burgers problem. This general formulation intro- 
duces the need for something like an entropy condition to achieve unique distribution 
solutions. For example, in the case of a negative delta-ring measure, we have a decreas- 
ing jump in m which introduces the possibility of a rarefaction solution going forward 
in time. In the classical weak solution formulation of Burgers equation, the entropy 
condition would select the rarefaction as the unique forward-time solution. Neverthe- 
less there exist other solutions, such as the outward-moving shock, that are bonafide 
distribution solutions, albeit ones that violate Lax's entropy condition whereby the 
speed of the shock should be faster than the characteristic speed ahead of it, and 
slower than the characteristic speed behind it. 

5. Conclusions. 

We have considered existence of radially symmetric, monotone decreasing solu- 
tions to the aggregation equation in the case of more singular potentials |x| a /a! for 
a in the range 2 — d < a < 1. We remind the reader that the problem with a > 1 
is known to be globally well-posed for measure data including the case without ra- 
dial symmetry and monotonicity |24j . For 2 — d < a < 2 we find that monotonicity 
is preserved, a feature that is not true for a > 2. Our results provide a rigorous 
framework for monotonicity behavior observed in numerical simulations of finite time 
blowup [44j [45] for radially symmetric data. The results also provide an understand- 
ing of the continuation of the solution after blowup. That understanding includes 
the result that one obtains instantaneous mass concentration for certain classes of L 1 
initial data including those observed as the asymptotic form of the blowup profile in 
numerical simulations (3H |3S] . The special case of the Newtonian potential results in 
a localization of the problem, reducing to a form of the inviscid Burgers equation on 
the half line. In particular for radially symmetric decreasing data, there is a unique 
classical solution of the Burgers problem for all time, resulting in a unique solution 
of the original density problem. This solution also retains its monotonicity. 

In contrast to the Newtonian potential, for the case 1 > a > 2 — d the ensuing 
velocity field is at best Holder continuous in time and our results arc less precise. For 

26 



example, uniqueness of solutions is still an open problem in this range, as is existence 
in the case of non-monotone, radially symmetric data. The existence problem is 
complicated by the fact that the velocity field is at best Holder continuous which makes 
it difficult to get convergence estimates for the flow map - something we use to prove 
existence of solutions in the case of monotone data. It is somewhat ironic that the more 
singular case of the Newtonian potential can be more easily solved - the velocity field 
is more singular, with a jump discontinuity. However the localization of the dynamics 
results in a better understanding of the problem. For the general nonlocal problem in 
the range 2 — d < a < 2, the monotonicity assumption allows for smoother estimates 
on the velocity field, namely Lipschitz estimates, which allow us to prove convergence 
of approximations and hence existence of a Lagrangian solution. In addition to the 
above problems for data with symmetry, the general problem of measure solutions with 
non-radially symmetric data is wide open. Some insights can be gained from recent 
work on special families of weak solutions. In the case of the Newtonian potential 
there exists a class of 'patch solutions' that are the time-dependent characteristic 
functions of of a domain in K d . These solutions have recently been observed [IT] to 
converge in finite time to a measure supported on a set of codimension one. Other 
works considers the analogue of vortex sheets for general aggregation equations with 
potentials that include both attraction and repulsion [751 GUI SI] ■ 

Acknowledgments. We thank the referees for many helpful comments. 

6. Appendix. 

6.1. Some general ODE results. 

In standard ODE textbooks such as [37J , it is proven that the flow map associated 
to a velocity field which is continuously diffcrentiable in both space and time is itself 
differcntiablc. In our case of interest the velocity field is continuously diffcrentiable in 
space but only continuous in time. We show here that the hypothesis of continuous 
differentiability in time can be replaced by a weaker assumption that holds true in 
our case. Only very minor modifications are needed compared to standard proofs 
found in ODE textbooks such as [27] . We will refer to [27] and we will indicate the 
necessary modifications to be made in the proof there. 

Let n C R d and J c (— oc, +oo) be two open sets. Suppose the function v : 
Six / satisfies the following: 

(HO) v is continuous on £1 x J. 

(HI) For every t 6 J, the function x *-> v(x,t) is continuously diffcrentiable on 17. 
(H2) Given compact sets f2i CC $7 and J\ CC J, there exists C > such that 

\Vv(x,t)\ < C 

for all (x,t) e!li x Ji. 
(H3) Given compact sets Oi CC SI and J\ CC J, and given e > 0, there exists 
S > such that 

\Vv(y,t) - Vv(x,t)\ < e 

for all x, y £ f2i satifying \x — y\ < S and for all t £ J\. 
In all the above Vi>(x,i) always stands for the derivative of v with respect to x. 

Theorem 6.1. Under the above hypothesis, given (xo,to) S Q X J there exists 
open sets f2o C f2 and Jq C J such that (xo,to) £ f^o x Jo an d a unique continuous 
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function a : il x J m> K d such that 

a(x,t) = x + I v(a(x, s), s)ds. (6-1) 



Moreover, given t G Jo, the mapping x i-> <r(a;,t) is continuously differentiable on flo 
and we have 

Vo-(x,t) = Id+ f Wv(o-(x,s),s)Vo-(x,s)ds. (6.2) 

J to 



Proof. The proof of Theorem 1.184, page 120 in [27] (or Theorem 1.261, page 
138 of the online version of |27j ) can be carried out with very minor modifications. 
Let us just mentioned where and how hypothesis (H3) is needed. In the proof of 
[2"7] the space X and Y are defined by X = C(b(t ,S) x B(xo 1 u/2),B(xq,v)) and 
Y = C b {b(t ,S) x B(x ,v/2) 1 L{R d ,R d )), where b(t ,S) and B(x ,v/2) denotes balls 
of radius 5 and vj2 and L(K d ,K rf ) denotes the set of linear transformations on R™. 
Both X and Y are endowed with the sup norm. Cb stands for continuous and bounded. 
The mapping # : X x Y — s- Y is defined by 

&((/), $)(x, t) = Id + [ Vv((f>(x,s),s)$(x,s)ds. 

In order to use the fiber contraction principle from [37J , we must verify that # is con- 
tinuous. From (H2) we easily obtain ||#(0, $i) - #(0, <& 2 )| < #(5 11*1 - $211 where 
A' = swp B ( Xo l/ /2)xb(t ,8) \^ v \- Hypothesis (H3) is needed in order to obtain continuity 
of # with respect to its first argument. To see this write 

#(<£i, $)(x, t) ~ #(0 2 , $)(x, i) = / (V«(0i (x, s), s) - Vv 2 (<j>(x, «),«)) s)ds. 

■/to 

(6.3) 

Using (H3) we see that |j#(0i, $) — #(02, < I ) )|| can be made as small as we want by 

choosing 0i and 02 close enough with respect to the sup-norm. 

□ 

Remark 4. Since Vi> is not assumed to be continuous with respect to time, 
the function t h-> V(r(x, t) is not necessarily continuously differentiable. However it 
is absolutely continuous as can been seen from (|6.2p . Therefore, given x G Qo, the 
function Y(t) — Va(x, t) is differentiable for almost every t G Jo and the differential 
equation 

Y'(t)=Wv{a(x,t),t) Y(t) 

holds almost everywhere in Jq. Then we can use Liouville Theorem (which is stated 
below) to deduce that 

— det S7a(x,t) = (div v)(a(x,t),t) det Va(x,t) (6.4) 
dt 

also holds almost everywhere in Jq . This of course implies that 



det S7a(x, t) = exp( I (div v)(a(x, s), s)ds) , 

J t 
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(6.5) 



which is the formula needed in our case (see (|2.30[) ). 

Theorem 6.2 (Liouville). Let A be dxd matrix and let Y(t) be a dxd time depen- 
dent matrix which is differentiable at t — to and satisfies Y'(to) = A Y(to). Then the 
function A(t) = det Y{t) is differentiable at to and satisfies A' (to) = (Tr A) A(to). 

Proof. See for example Hartman [40] , □ 

6.2. Some general Lemmas. 

A proof of the following Lemma can be found in [551 Lemma 3.11]. 

Lemma 6.3. Let T,S :R d -^R d be two Borel maps. Also take p £ V 2 (R d ). Then 

W 2 (S#p,T#p) < \\S-T\\ L 

Lemma 6.4. Suppose that p 6 V(R d ) has compact support and suppose that 
Vrd^ 1 ) 3 p n converges narrowly to p. Then p also belongs to VnD(^- d )- 

Proof. Let R be a rotation of R d and let / e C(R d ). Then / o R g C(R d ) and 
/ foRdpn = J fdp n . Taking limits we see that J foRdp = J fdp so that p £ A4n(R d ). 
To prove p is decreasing fix < r\ < r 2 and take disjoint small rings Aj = {r.j — r)j < 
\x\ < rj + Sj},j = 1,2 having the same volume. We may assume p(dAj) = 0. 
Then there exist continuous functions f\ > xa 1 and f 2 < Xa 2 with disjoint supports 
such that \p(Aj) — J fjdp\ < e. By hypothesis we have inf f\dp n > J f 2 dp n , so that 
p(A 2 ) + 2e < p(Ax). Then shrinking e, A x and A 2 shows that p e V RD (R d ). □ 

Lemma 6.5. Let a > 2 — d and suppose y : [0,+oo) — > [0,+oo) is an absolutely 
continuous function satisfying — y(t) Q_1 < y'(t) < for almost every t £ [0, +00) for 
which y{t) > 0. If a < 1 then y(t) is Holder continuous. To be more precise: 

_((2- Q )(t-s))^ <y(t)-y(s)<0 
for allO < s <t. If a > 1 then y(t) is Lipschitz continuous. To be more precise: 

-y(0) a - 1 (t~ S )<y(t)-y(s)<0 

for allO < s <t. 

Proof. The case a > 1 is trivial since the inequality ~J/(^) Q_1 < y'{t) < together 
with the non-negativity of y implies — 2/(0) Q_1 < y'{t) < 0. We now prove the Lemma 
for 2 — d < a < 1. For almost every t > for which y(t) > we have —(2 — a) < 

ft (y(*) 2_Q ) < °- lt is then clcar that ~( 2 - «)(* - s ) < y(t) 2 ~ a - y(s) 2 ~ a < for 

all < s < t. But because of the convexity of the function r M> r 2 ~ a we have that 
(y(s) — y{t)) 2 ~ a < y(s) 2 ~ a — y(t) 2 ~ a , which gives the result.Q 
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